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ABSTRACT 

1. Detecting aquatic macroorganisms with environmental DNA (eDNA) is a new survey method 
with broad applicability. However, the origin, state, and fate of aqueous macrobial eDNA - 
which collectively determine how well eDNA can serve as a proxy for directly observing 
organisms and how eDNA should be captured, purified, and assayed - are poorly understood. 

2. The size of aquatic particles provides clues about their origin, state, and fate. We used 
sequential filtration size fractionation to measure, for the first time, the particle size distribution 
(PSD) of macrobial eDNA, specifically Common Carp (hereafter referred to as Carp) eDNA. We 
compared it to the PSDs of total eDNA (from all organisms) and suspended particle matter 
(SPM). We quantified Carp mitochondrial eDNA using a custom qPCR assay, total eDNA with 
fluorometry, and SPM with gravimetric analysis. 

3. In a lake and a pond, we found Carp eDNA in particles from >180 to <0.2 um, but it was most 
abundant from 1-10 um. Total eDNA was most abundant below 0.2 um and SPM was most 
abundant above 100 um. SPM was <0.1% total eDNA, and total eDNA was <0.0004% Carp 
eDNA. 0.2 um filtration maximized Carp eDNA capture (85 %±6%) while minimizing total (i.e., 
non-target) eDNA capture (48%±3%), but filter clogging limited this pore size to a volume <250 
mL. To mitigate this limitation we estimated a continuous PSD model for Carp eDNA and 
derived an equation for calculating isoclines of pore size and water volume that yield equivalent 
amounts of Carp eDNA. 

4. Our results suggest that aqueous macrobial eDNA predominantly exists inside mitochondria or 
cells, and that settling plays an important role in its fate. For optimal eDNA capture, we 
recommend 0.2 um filtration or a combination of larger pore size and water volume that exceeds 
the 0.2 um isocline. In situ filtration of large volumes could maximize detection probability 
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when surveying large habitats for rare organisms. Our method for eDNA particle size analysis 
enables future research to compare the PSDs of eDNA from other organisms and environments, 
and to easily apply them for ecological monitoring. 

KEY-WORDS: environmental DNA, particle size analysis, sampling methods, ecological 
monitoring, genetic monitoring, aquatic ecosystems, rare species 

INTRODUCTION 

Environmental DNA (eDNA) is DNA extracted from bulk environmental samples (e.g., soil, 
water, air) without isolating target organisms or their parts from the sample. The concept and the 
term both originate from microbiology (Ogram et al. 1987) where the target DNA in 
environmental samples is from abundant live and dead microbes. In contrast, macrobial eDNA is 
the DNA of large organisms such as animals or plants that occurs in environmental samples. 
Although macrobial eDNA has been studied since 1991 in fields such as human forensics 
(Hochmeister et al. 1991), agricultural transgenics (Widmer et al. 1997), paleogenetics 
(Willerslev et al. 2003), and fecal pollution source tracking (Martellini et al. 2005), it was only 
in 2008 that eDNA was first used to detect aquatic macrofauna (Ficetola et al. 2008). Aqueous 
macrobial eDNA has garnered particular interest as a simple and sensitive alternative to directly 
observing rare aquatic macrofauna. Despite a recent burst of papers demonstrating its use for 
detection of diverse macrofauna, we have much to learn about the origin (e.g., primary 
physiological source), state (e.g., intra- or extra-cellular), and fate (e.g., suspension time) of 
aqueous macrobial eDNA (Lodge et al. 2012). These three domains collectively determine how 
well eDNA can serve as a proxy for directly observing organisms, and how eDNA should be 
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captured to make robust inferences about organism presence or abundance. In this paper we 
determine the size distribution of particles containing eDNA and draw size-based inferences 
about the state and fate of eDNA. 

Martellini et al. (2005) first demonstrated that natural waters contain mitochondrial DNA 
(mtDNA) from vertebrates inhabiting the watershed. Their research focused on terrestrial 
mammals for which aqueous eDNA was presumably of allochthonous fecal origin. 
Autochthonous aqueous eDNA from aquatic macrofauna was subsequently described for 
amphibians (Ficetola et al. 2008), fish (Jerde et al. 201 1), mammals, crustaceans, insects 
(Thomsen et al. 2012a), birds (Thomsen et al. 2012b), and reptiles (Piaggio et al. 2013) in lentic, 
lotic, and marine waters. eDNA-based surveys for rare fish or amphibians are more sensitive 
than traditional methods (Jerde et al. 201 1; Thomsen et al. 2012b; Dejean et al. 2012) and eDNA 
disappears from surface water within 1 to 25 days of organism absence (Dejean et al. 201 1; 
Thomsen et al. 2012a; b). Aqueous macrobial eDNA concentration is very low (up to 33 mtDNA 
copies mL" 1 ) and correlates with organism density, at least for amphibians in small ponds and 
streams (Thomsen et al. 2012a; Pilliod et al. 2013) and fish in experimental ponds (Takahara et 
al. 2012). Aqueous macrobial eDNA has been captured (i.e., concentrated from water) by 
precipitation/centrifugation (Martellini et al. 2005; Ficetola et al. 2008), chromatography 
(Douville et al. 2007), lyophilization (Pote et al. 2009), ultrafiltration (pore size <0.1 um; 
Takahara et al. 2012), and microfiltration (pore size >0.1 um; Martellini et al. 2005; Goldberg et 
al. 201 1; Jerde et al. 201 1; Minamoto et al. 201 1). 
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Despite the considerable variation in methods used to capture aqueous macrobial eDNA (Table 
SI), we are unaware of any systematic comparison between methods or description of the 
captured particles. The size of aquatic particles helps determine their characteristics and 
interactions with organisms, other particles, and the environment (Burd & Jackson 2009). Thus 
the size of aquatic particles containing macrobial eDNA is foundational to our understanding of 
the origin, state, and fate of this material. For example, large particles sink faster than small 
particles (Maggi 2013), so eDNA-based surveys aimed at determining very recent and local 
organism presence might need to target the upper end of the eDNA-containing particle size 
distribution (PSD). The PSD of aqueous macrobial eDNA can also inform decisions about 
sampling effort (e.g., water volume) and particle capture size (e.g., filter pore size). Currently 
these decisions are made using only precedent, trial and error, and logistical constraints. Finally, 
the PSD of macrobial eDNA may help mitigate one of the most troublesome issues associated 
with genetic testing for low copy number DNA in environmental samples: sample interference 
(Schrader et ah 2012). Non-target material suspended in water, including non-target DNA 
(Thompson et ah 2006), can interfere with eDNA recovery (Lloyd et ah 2009) and genetic 
assays (Hedman & Radstrom 2013). If the PSD of macrobial eDNA differs from that of non- 
target material then size-based, selective enrichment of target eDNA may be possible. 

Here we produce the first measurement of PSD for aqueous macrobial eDNA. We sequentially 
filtered water from a lake and pond inhabited by Common Carp {Cyprinus carpio Linnaeus, 
1758, hereafter Carp) to create a series of size fractions containing the particles retained at each 
size, from >1 80 um to <0.2 um. In aquatic microbiology, 0.2 um is the filter pore size below 
which aqueous eDNA exists as extracellular molecules free in solution (Matsui et ah 2004; 
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Maruyama et al. 2008), and the vertebrate mitochondrion measures 0.2-8 um (Flindt 2006 p. 
254). Thus our size fractions spanned from extracellular and extraorganellar DNA (hereafter 
extramembranous; <0.2 um) to particles larger than most vertebrate cells (>1 80 um). For each 
size fraction, we measured the concentration of Carp eDNA, total eDNA (DNA of any type), and 
suspended particle matter (SPM). This allowed us to compare PSDs among these three nested 
types of aquatic matter and evaluate the opportunity for size-based target enrichment. 

MATERIALS & METHODS 

Target species 

We selected Carp as the target macrobial species because fish species, including Carp, have been 
the target of many previous studies of aqueous macrobial eDNA (Dejean et al. 201 1; Jerde et al. 
201 1; Minamoto et al. 201 1; Takahara et al. 2012, 2013; Thomsen et al. 2012a; b; Collins et al. 
2012; Wilcox et al. 2013). In addition, Carp is listed at number 30 on the IUCN list of the 
world's worst invasive species (Lowe et al. 2004). Because Carp inhabit many continents, our 
method for particle size analysis can be applied globally to assess the consistency of one species' 
eDNA PSD across environmental conditions. 

Species-specific qPCR assay for eDNA quantification 

We designed a species-specific SYBR Green qPCR assay to measure the concentration of only 
Carp eDNA in environmental samples. The qPCR primers were manually designed in an 
alignment of all available mitochondrial control region (D-loop) sequences from NCBI's 
GenBank for Carp and the most closely related species that potentially co-occur in North 
America: Goldfish (Carassius auratus), Grass Carp (Ctenopharyngodon idella), Black Carp 
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{Mylopharyngodon piceus), Crucian Carp (Carassius carassius), and Bigheaded Carp 
{Hypophthalmichthys spp.). We conducted further in silico testing of species specificity using 
NCBI Primer-BLAST tool. A set of candidate primers was selected for in vitro testing on tissue- 
derived total genomic DNA from Carp, Goldfish, Grass Carp, Black Carp, and Bigheaded Carp. 
We selected one primer set (CcrpDlpF: GAGTGCAGGCTCAAATGTTAAA; CcrpDlpR: 
GTAAGGATAAGTTGAACTAGAGACAG; 146 bp amplicon) that demonstrated efficient 
amplification from Carp and no amplification from any of the non-target species. Finally, we 
conducted in situ testing of this assay using eDNA samples from natural water bodies with and 
without Carp. The assay did not amplify eDNA from a small lake without Carp (Knop Lake, 
Indiana, USA). The assay amplified eDNA from two large lakes with Carp (Lakes Michigan and 
Huron, Michigan, USA), and Sanger sequences (ABI 3730x1, Applied Biosystems) of these 
qPCR products matched the target amplicon from Carp. The DNA melting temperature (T m ) of 
all qPCR products from all eDNA samples matched that from tissue-derived Carp DNA, as 
determined by a melting curve. 

During development and application of the assay, we used 20 uL reactions containing 10 uL of 
Power SYBR Green Master Mix (Life Technologies), 4 uL of DNA, 2 uL of 4 ug*uL~' non- 
acetylated BSA (Life Technologies), and 2 uL of each 3 uM primer. We performed all reactions 
on an Eppendorf Mastercycler ep realplex2 S thermocycler with the following reaction 
conditions: 95°C for 10 min.; 55 2-step cycles of 95°C for 15 sec. and 60°C for 1 min., 95°C for 
15 sec, 60°C for 15 sec, and a 20 min. melting curve up to 95°C for 15 sec. Fluorescence data 
collection occurred during the 1 min. 60°C step and during the melting curve. 
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We used absolute quantification qPCR with a standard curve to quantify Carp eDNA. We made 
the standard curve DNA by amplifying the entire mtDNA control region (1022 bp) (Liu & Chen 
2003) from tissue-derived Carp DNA. After confirming a single gel electrophoresis band 
approximately 1000 bp in size we purified this PCR amplicon using ExoSAP and quantified it 
using 5 uL of PCR product with a Qubit fluorometer (Invitrogen) and the Qubit dsDNA High 
Sensitivity kit. We converted from DNA weight to DNA copies using the median molecular 
weight of the 95% consensus 1022 bp amplicon sequence from all C. carpio mitogenomes on 
GenBank (631214 g*mole -1 ) as calculated by OligoCalc (Kibbe 2007). We stored single use 
7.7E+06 copies*uL _1 aliquots of this standard DNA at -20°C. qPCR standard curves were 
prepared by serial 10-fold dilution from 7.7E+06 copies*uL _1 down to 1 copy*uL _1 in low TE 
buffer (lOmM Tris, O.lmM EDTA) and stored at 4°C while in use (Bellete et al. 2003). All 
manipulations of standard DNA, eDNA, and reagents for qPCR used low bind tubes and low 
bind aerosol barrier pipette tips (Ellison et al. 2006). 

Water bodies 

We collected water from two lentic water bodies, selected to provide a strong contrast in 
simplicity of the physical environment, in macrofaunal diversity, and in abundance of the target 
species (Carp). At Potawatomi Zoo in South Bend, Indiana, USA (41.670629, -86.216634), we 
sampled a small (0.06 ha, maximum depth 2 m) concrete-lined outdoor pond with flow-through 
treated municipal water that contained approximately 500 large adult Carp. We expected that this 
extremely high Carp density would provide sufficiently high carp eDNA concentration to allow 
quantification in every size fraction following sequential filtration. The more natural, lower Carp 
density environment was St. Mary's Lake at the University of Notre Dame, Indiana, USA 
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(41.701497, -86.244312), a natural spring-fed mesotrophic lake (10.7 ha, maximum depth of 9.1 
m). We lack a quantitative estimate of Carp density for this lake but based on our many visual 
observations, the density of Carp is far lower in the lake than in the zoo pond. 

Water collection and sampling 

In December 201 1 ("winter") we conducted our primary sampling event for each water body. In 
April 2012 ("spring") we conducted a supplemental sampling event at St. Mary's Lake because 
Carp eDNA concentrations measured in winter were very low (see Results and Discussion). For 
each sampling event we pumped 12 L of water from <1 m depth into an autoclaved carboy using 
autoclaved silicone peristaltic tubing with internal diameter of 9.5 mm (Masterflex L/S 36) and a 
portable, battery-powered peristaltic pump (Alexis, Pegasus Pump Company, Bradenton, Florida, 
USA). At the pond we pumped water from the east end of the pond where Carp were aggregated 
by extending the tubing approximately 2 m from shore. At the lake we pumped water from a 
small embayment on the south end as we rowed a bleach-decontaminated boat to and from a 
point approximately 100 m from shore. Carboys were transported to the laboratory within 
approximately 20 minutes and water was distributed into autoclaved 300 mL bottles. We stored 
the 300 mL bottles on ice until each was used for particle size fractionation, which began 
immediately and finished within 10 hours. 

Particle size fractionation 

Filtration can bias particle size analysis by retaining particles smaller than the filter pore size or 
passing particles larger than the pore size (Droppo 2006). However, our analysis required 
chemical assay of the SPM to quantify its subcomponents, total eDNA and Carp eDNA. Thus 
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particle sizing methods that do not require fractionation (e.g., electrical, optical, or flow 
cytometry particle counting) were not applicable (Reynolds et al. 2010). Filtration was also the 
most appropriate method given our aim to use PSDs for informing optimal eDNA capture - 
which is typically done by filtration. For these reasons we used sequential filtration size 
fractionation (Figure 1) to generate the PSDs for SPM, total eDNA, and Carp eDNA. We 
followed the recommendations of Droppo (2006) to minimize filtration artifacts: low vacuum 
pressure, small water volume, and filters with precise pores. We used nylon net filters (180, 100, 
and 60 um mesh size; Millipore) and track-etched polycarbonate (PCTE) filters (20, 10, 1, 0.2 
urn pore size; GE Osmonics) because these provide the sharpest separation between particle 
sizes and the best agreement between nominal and effective pore size (Sheldon 1972; Droppo 
2006). By contrast, 'tortuous path' type filters such as glass-fiber and cellulose membrane 
actually have a wide array of effective pore sizes in every filter, rather than well-defined pores 
(see Droppo 2006 for a detailed discussion). Thus our PSDs, while accurate, may not translate 
directly to the nominal pore sizes of tortuous path filters, and we recommend the use of PCTE 
filters for future particle size analysis of aqueous macrobial eDNA. 

From each of the ten (winter) or three (spring) bottles, 250 mL was sequentially filtered through 
47 mm diameter filters of decreasing pore size using gentle vacuum pressure (<64 cm Hg). We 
used 47 mm diameter magnetic filter funnels (Pall) and 1 L polypropylene vacuum flasks 
(Nalgene) that were decontaminated between uses by immersion in a 10% bleach solution 
followed by copious rinsing with reverse osmosis treated water. To avoid inaccurate size 
fractionation caused by filter clogging (which we observed in a pilot experiment), we used only 
150 mL of the 1 um filtrate for filtration through the 0.2 um filter. 15 mL of the 0.2 um filtrate 
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was transferred to a 50 mL centrifuge tube containing 33.5 mL of 100% ethanol and 1.5 mL of 3 
M sodium acetate then frozen at -20°C for subsequent eDNA capture following the precipitation 
procedure of Ficetola et al. (2008) with a centrifuge force of 3220 g. We placed used filters in 
CTAB buffer (Coyne et al. 2005) and stored them at -20°C until subsequent DNA extraction. 

For the primary sampling event (winter) only, we conducted size fractionation for gravimetric 
analysis at the same time and in the same manner as for DNA analysis except all filters were 
individually pre-desiccated, pre-weighed and stored in individually labeled petri dishes. 
Following filtration, a filter disc was carefully returned to its petri dish retentate side up and 
placed in a 60°C desiccation chamber for 24 hours prior to gravimetric analysis. The final size 
fraction (<0.2 um) was not included in gravimetric analysis because we were unable to weigh the 
tiny pellet produced by precipitation/centrifugation of the 0.2 um filtrate. 

Assaying eDNA: contamination precautions 

We conducted eDNA extraction, total eDNA quantification, and qPCR setup in a strictly pre- 
PCR laboratory located on a separate floor and opposite end of the building from our post-PCR 
laboratory. eDNA extraction and total eDNA quantification were performed on a bench used 
exclusively for those purposes and qPCR setup was performed on a different bench used 
exclusively for that purpose which is located on the opposite end of the pre-PCR laboratory 
(Mifflin 2007). A DNA extraction negative control was included for every filter pore size and for 
every batch of extractions. Each qPCR plate included three no template control (NTC) reactions. 

Assaying eDNA: eDNA extraction 
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We extracted eDNA from filters and precipitation pellets following the CTAB extraction of 
Coyne et al. (2005). PCTE filters dissolve during the extraction and nylon net filters were 
shredded in their extraction tubes using flame- and bleach-sterilized scissors. Final eDNA pellets 
were re-suspended in 100 uL of low TE buffer. 

Assaying eDNA: total eDNA quantification (Qubit) 

We quantified total eDNA using 5 uL of eDNA extract in the Qubit dsDNA High Sensitivity kit 
and a Qubit fluorometer (Invitrogen). The limit of detection for this method is 0.5 ng*mL _1 in the 
200 uL assay tube (corresponding to 0.02 ng*uL~' in our eDNA extracts) thus we set the eDNA 
concentration at 0.02 ng*uL~' for eDNA extracts with concentrations too low for Qubit 
measurement. This reflects our assumption that eDNA existed at every size fraction, even if the 
sample volume limitations of our size fractionation produced some extract concentrations too 
low for Qubit detection. 

Assaying eDNA: Carp eDNA quantification (qPCR) 

We measured Carp eDNA concentration in each eDNA extract using triplicate reactions of the 
qPCR assay performed as described above. To minimize variation between replicate reactions 
caused by imperfect pipetting of small eDNA extract volumes, we combined eDNA extract and 
master mix (all other reagents) for three reactions into one tube then dispensed to three plate 
wells using an electronic repeating pipette (Xplorer 5-100 uL, Eppendorf). Each qPCR plate 
included a five point standard curve from 3.1E+04 copies*reaction~' down to 3 copies*reaction _1 . 
The fluorescence threshold for each plate was automatically determined by the Eppendorf 
realplex software using the default Noiseband setting. The fluorescence baseline was calculated 
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for every reaction individually using the default Automatic Baseline setting of the Eppendorf 
realplex software. Every amplification profile and melting curve profile was visually examined 
to confirm exponential amplification and a T m matching that of standard curve reactions. 
Following the recommendation of Ellison et al. (2006) for low copy number qPCR, we 
calculated concentrations for each reaction, assigning zero concentration to non-detect reactions 
and averaging concentration across the three technical replicates for an eDNA extract. In some 
reactions the measured reaction copy number was slightly below one so we rounded all reaction 
copy numbers up to the next largest integer. 

Assaying SPM (gravimetric analysis) 

Desiccated filters containing captured particles were weighed using a microbalance. To eliminate 
the effect of static charge buildup on microbalance readings, we passed every filter through two 
2U500 polonium-210 ionizers held in a Staticmaster positioner (Amstat Industries, Inc., 
Glenview, Illinois, USA). 

Data conversion, analysis, and modeling 

To compare the PSD shape of Carp eDNA, total eDNA, and SPM, we converted concentration in 
each size fraction to percent of total (Figure 2). This assumed that the final fractionation step 
(precipitation for Carp eDNA and eDNA, 0.2 urn filtration for SPM) captured all remaining 
matter for each of the three types, respectively. To evaluate the optimal capture of Carp eDNA, 
we converted concentration of Carp eDNA, total eDNA, and SPM to cumulative percent 
captured (capture efficiency, E) at each fractionation step. To allow the most general evaluation 
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possible from our experiment, we pooled the lake and pond capture efficiency data (winter only) 
without making any assumptions about the similarity or difference between the two water bodies. 
To estimate Carp eDNA capture efficiency as a continuous function of filter pore size (i.e., a 
continuous PSD for Carp eDNA) we modeled our discrete capture efficiency data with the 
complementary cumulative distribution function (CCDF) of the Weibull distribution, a flexible 
distribution that was originally developed and is widely used to describe PSDs (Rosin & 
Rammler 1933; Zobeck et al. 1999; see Appendix SI for details). Finally, we used this 
continuous model to calculate 'isoclines' of filter pore size and filtered water volume where the 
Carp eDNA PSD estimates identical captured amounts of Carp eDNA. 

To calculate the proportion of total eDNA that was made up by Carp mtDNA, we converted the 
qPCR-measured Carp eDNA copy number to weight of Carp mitogenomes. This was done using 
the molecular weight of one arbitrarily chosen Carp mitogenome from NCBI GenBank 
(10244181 g^mole" 1 ; accession number: JN105352.1) as calculated by OligoCalc (Kibbe 2007). 

RESULTS 

All extraction negative controls and qPCR NTCs tested negative for Carp eDNA. In both water 
bodies, particles containing Carp eDNA were detected in all of the size fractions, from > 180 um 
to < 0.2 um, and Carp eDNA was most abundant in the 1-10 um particle size fraction (Figure 2). 
On average there was more Carp eDNA in particles larger than 0.2 um for both water bodies 
(pond: 10 of 10 replicates; lake: 7 of 10 in winter, 3 of 3 in spring) although this difference was 
significant only for the pond (Chi-square goodness of fit test; pond: % (1, n=10) = 10.00, p < .01; 

lake-winter: % 2 (1, n=10) = 1.60, p = 0.21; lake-all: % 2 (1, n=13) = 3.77, p = 0.05). The average 
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percent of Carp eDNA in particles larger than 0.2 urn was 93%±1% (mean ± s.e.m; pond), 
75%±13% (lake-winter), and 78%±10% (lake-all). These results demonstrate that Carp eDNA 
can exist in very large particles and suggest that extramembranous Carp eDNA (particles <0.2 
um) represents a small portion of the total size distribution of particles containing Carp eDNA. 
In contrast, total eDNA was most abundant in the <0.2 um size fraction, and SPM was most 
abundant in the 100-180 um size fraction (pond) or the >1 80 um size fraction (lake) (Figure 2). 

The average percent of SPM that was total eDNA was 0.01%±0.001% (mean ± s.e.m.; pond) and 
0.1%±0.01% (lake) when summed across all size fractions. The average percent of total eDNA 
that was Carp mitochondrial eDNA was 0.0004%±0.0001% (pond) and 
0.0000004%±0.0000001% (lake) when summed across all size fractions. These results 
demonstrate that total eDNA and SPM were nearly 100% non-target material with respect to the 
target, Carp mtDNA. 

Filtration with a 0.2 um pore size maximized Carp eDNA capture efficiency (85%±6%; 
cumulative data pooled across lake and pond) while minimizing total eDNA capture efficiency 
(48%±3%). Unfortunately, 0.2 um pore size filters clog with very small throughput volumes 
(e.g., we were unable to filter 250 mL of 1.0 um filtrate through our 0.2 um filter). Maximum 
throughput volume can be increased by either increasing filter surface area or increasing filter 
pore size. Our data cannot provide guidance for increasing filter surface area, but we used the 
continuous PSD of Carp eDNA to estimate the combinations of pore size and volume that yield 
identical amounts of Carp eDNA. The estimated scale and shape parameters (X and k) from the 
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Weibull CCDF model (i.e., the continuous PSD) of Carp eDNA were 1 1.85 and 0.351 and the 
estimated Carp eDNA capture efficiency at 0.2 um (E0.2) was 78.8% (Figure SI). 



Thus the equation for estimating isoclines relative to 0.2 um filtration is derived as follows: C* 
V0.2 *Eo.2 = C*V X *E X (eqn 1), where C is the environmental concentration of the target, V is 
throughput volume, and E is target capture efficiency. Solving for V x yields V x = (Vo.2*E 0 .2) / E x 
(eqn 2). V0.2 is whatever throughput volume the researcher determines is feasible for a 0.2 um 
pore size filter, E0.2 is target capture efficiency at 0.2 um, and E x is the continuous PSD of the 
target, yielding the final equation: 

V x = — tttt (eqn 3) 

Figure 3 plots isoclines using the continuous PSD for Carp eDNA and several hypothetical 
values of Vo.2- 



DISCUSSION 

Methodological considerations 

As this paper represents the first investigation of particle size for aqueous macrobial eDNA, we 
first consider some methodological points that will be important should other scientists use our 
approach. The per reaction copy numbers of Carp eDNA measured by qPCR were very low, with 
a heavily right-skewed distribution and a mode of zero for both water bodies (Figure S2). This 
result is consistent with detailed empirical studies of qPCR applied to low copy number DNA 
(i.e., < 100 target copies per reaction) (Ellison et al. 2006). Research using qPCR to quantify 
eDNA should report how 'non-detects' (negative reactions from technical replicate sets where at 
least one reaction was positive) are handled. As demonstrated by Ellison et al. (2006), these non- 
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detect reactions should be assigned a zero target concentration and included when averaging 
across technical replicates, as this provides better estimates of unknown sample concentration. 

Even with this low copy number protocol applied, the frequency of negative reactions (79%) and 
non-detects (22%) in the lake-winter sampling event was high. The large sampling error (Figure 
2) caused by these extremely low target DNA levels led us to conduct the supplemental spring 
sampling at the lake, because we expected higher fish activity in the spring would increase Carp 
eDNA concentration. Indeed, the frequency of negative reactions and non-detects dropped to 
44% and 19%, respectively, in the lake-spring sampling event. The spring PSD from the lake 
supported our conclusions that Carp eDNA exists predominantly in particles larger than 0.2 um 
and is most abundant between 1 and 10 um (Figure 2). The supplemental spring sampling event 
also confirmed that total Carp eDNA concentration in the lake (summed across all size fractions) 
had increased four-fold since winter (from 1726±602 to 7144±1 165 copies*L _1 ). These data 
represent the first description of seasonal change in eDNA concentration for aquatic macrofauna 
and underline the importance of reporting eDNA sampling dates, especially when comparing 
eDNA sampling with organism sampling (Thomsen et al. 2012a). 

State, and fate of aqueous macrobial eDNA 

The majority of Carp eDNA was found in particles larger than 0.2 um in both the pond (93%) 
and the lake (75%-78%), and Carp eDNA was most abundant in the 1-10 um particle size 
fraction (Figure 2). In contrast, total eDNA was most abundant in the <0.2 um size fraction and 
SPM was most abundant in the 100-180 um (pond) or the <180 um size fraction (lake) (Figure 
2). Future research from other seasons, environments, and species will determine the generality 
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of these results, but even at this early stage it seems appropriate to consider how the PSD can 
shed light on the state and fate of aqueous macrobial eDNA. 

Aqueous eDNA from living macrofauna most likely originates as urine and feces, epidermal 
tissues and secretions, and reproductive cells and fluids. Much of this source material (e.g., 
feces) will enter the water column as large particles (>1000 urn), yet we generally found little 
Carp eDNA in particles larger than 60 um. Therefore we conclude that these sources of eDNA 
must rapidly settle or break apart. Indeed, particle size is a major determinant of settling velocity, 
even for complex biological and biomineral aggregates (Maggi 2013). It is also well-documented 
that animal feces contain viable epithelial cells (10 _1 -10 6 cells*g _1 ) and large amounts of the 
animal's mtDNA (10 _1 -10 7 copies*g _1 ) (reviewed in Caldwell et al. 201 1) and that feces from 
aquatic macrofauna rapidly sink (Robison et al. 1981; Wotton & Malmqvist 2001). Thus at least 
one obvious source of aquatic macrobial eDNA appears to spend very little time suspended in 
the water column. We measured PSDs in still water whereas particle settling should be slower in 
flowing water, such as streams. 

On the small end of the size spectrum there was generally little extramembranous Carp eDNA 
(<0.2 um), at least as indicated by the 146 bp mtDNA region our qPCR assay targeted. One 
might predict an accumulation of smaller particles as discharged materials decompose, ultimately 
leading to an abundance of free eDNA. However, free eDNA may be exceptionally susceptible 
to hydrolysis given its full exposure to microbial extracellular enzymes that are abundant in 
aquatic systems (Matsui et al. 2001). Thus the abundance of Carp eDNA in the 1-10 um fraction 
likely reflects a state of eDNA persisting within small cells and/or mitochondria. We cannot 
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exclude the possibility that Carp eDNA in particles larger than 0.2 um is actually 
extramembranous DNA aggregated or adsorbed onto larger particles (Burd & Jackson 2009; 
Suzuki et al. 2009), but the fact that animal mitochondria range in diameter from 0.2-1.2 um and 
in length from 1-8 um (Flindt 2006) suggest mitochondria as a parsimonious explanation for the 
pattern we observed. Indeed, during the regular apoptotic shedding of epithelial cells, intact 
mitochondria are released within apoptotic bodies and mtDNA is protected from the 
endonuclease degradation that rapidly degrades nuclear DNA (Murgia et al. 1992; Tepper & 
Studzinski 1993). For aquatic animals, this process will release whole mitochondria into the 
water column where their double membrane may resist lysis while the mitochondrial nucleoid 
further protects mtDNA (Rickwood & Chambers 1981). 

The PSD of Carp eDNA indicates that settling, in addition to degradation, plays an important 
role in the fate of aqueous macrobial eDNA. Particles larger than 1 um settle in natural waters 
(Isao et al. 1990) and 71%±5% of particles containing Carp eDNA exceeded this size (Figure 2). 
SPM in the 1-100 um size range settles at terminal velocities between 0.4 mm^hr" 1 and 40 m*hr" 
, depending on size, shape, and composition (Maggi 2013). Thus it appears that most aqueous 
macrobial eDNA is in a constant state of downward flux. Most research on the fate of aqueous 
macrobial eDNA has not considered settling (Martellini et al. 2005; Kortbaoui et al. 2009; 
Dejean et al. 201 1; Thomsen et al. 2012a), except one experiment that applied water circulation 
to eliminate it (Thomsen et al. 2012b). Our findings emphasize that degradation of DNA 
molecules and particle settling combine to reduce aqueous eDNA concentration over time and 
future research should try to measure both processes. Past eDNA degradation experiments may 



19 



Downloaded from http://biorxiv.org/on September 18, 2014 



have underestimated the persistence of aqueous eDNA given that settled particles can return to 
the water column through bioturbation and water turbulence (Bloesch 1995). 

We did detect a non-trivial amount of Carp eDNA in particles smaller than 1 um (Figure 2). This 
fraction of aqueous macrobial eDNA is likely to remain in suspension indefinitely until the DNA 
molecules degrade (Isao et al. 1990). Water could advect eDNA in this state over long distances, 
especially in lotic systems, strong currents, and when free eDNA is protected via adsorption to 
other submicron particles (Saeki et al. 201 1). For example, transgenes from genetically modified 
corn were detected in river water up to 82 km downstream of a corn cultivation plot (Douville et 
al. 2007), and plant DNA fragments as long as 1000 bp are detectable in groundwater (Pote et al. 
2009). The stability and advection of free DNA in aquatic environments is sufficient that 
scientists use synthetic DNA as a hydrologic tracer (Sabir et al. 2002; Foppen et al. 201 1). 
Future work examining the transport and persistence of naturally occurring aqueous macrobial 
eDNA across long distances will improve our understanding of the spatiotemporal window for 
inferring organism presence. The size of this window may depend on the particle size fraction 
that is tested for eDNA. 

Optimal capture of aqueous macrobial eDNA 

Filtration with a 0.2 um pore size was the best strategy for maximizing Carp eDNA and 
minimizing non-target eDNA (Figure 2). However, the continuous PSD of Carp eDNA that we 
estimated provides simple guidance for overcoming the volume limitations of 0.2 um filtration 
by increasing pore size and water volume along an isocline (Figure 3). These isoclines also 
highlight the potential to exceed the yield of 0.2 um filtration by using large pore sizes (e.g., >50 
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urn) and filtering volumes above the isocline. Given the basic processes discussed above 
regarding the origin, state, and fate of macrobial eDNA, we suspect that these eDNA capture 
recommendations are applicable beyond Carp. However, as researchers generate eDNA PSDs for 
additional species and environments, they can easily be compared via Weibull CCDF modeling 
and our simple equation for calculating isoclines. 

The heterogeneous distribution of macrobial eDNA among replicate water samples observed in 
this study (Figure 2) and others (Jerde et al. 201 1; Dejean et al. 2012; Thomsen et al. 2012a; 
Pilliod et al. 2013), suggests an extensive search mode (i.e., large pore size and large water 
volume) as the default strategy when conducting eDNA-based monitoring of rare species in large 
habitats. Devices for in situ filtration of large volumes, such as inline filters or plankton nets, 
could maximize water volume and habitat coverage. Our data indicate that capture methods with 
extremely small size cutoffs (e.g., precipitation, centrifugation, ultrafiltration) are generally 
unnecessary compared with microfiltration (Figure 2). The >50% increase in non-target eDNA 
capture (Figure 2) and the severe volume limitations of these methods further reduce their value 
relative to microfiltration. 

Because macrobial eDNA is only a tiny fraction of the eDNA and SPM in water, the risk of 
interference from non-target eDNA and non-DNA inhibitors increases with sample volume 
(McDevitt et al. 2007; Hata et al. 201 1; but see Gibson et al. 2012). Filtration reduces the risk 
from non-target eDNA by letting 50% or more pass through (Figure 2), and pre-filtration could 
reduce the risk from non-DNA inhibitors as SPM was most abundant on 100 or 180 urn filters 
(Figure 2). In contrast to non-target eDNA, the interference risk from SPM is mitigated by DNA 
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extraction and purification, thus we recommend pre-filtration only if it substantially increases the 
throughput volume of one's capture filter. For example, Cary et al. (2006, 2007) applied this 
'band-stop' filtration approach for eDNA-based surveillance of the invasive diatom 
Didymosphenia geminate. 

There is broad interest in research on the reliability of eDNA-based ecological monitoring 
(Sutherland et al. 2013; Zhan et al. 2013; Schmidt et al. 2013). Particle size analysis, such as that 
presented here, is one of the many important research avenues needed to better understand the 
origin, state, and fate of aqueous macrobial eDNA. Our first description of particle size for 
aqueous macrobial eDNA provides immediate guidance for practitioners and a tested method for 
researchers. We look forward to future findings that will undoubtedly expand on this initial 
description by investigating the organismal and environmental determinants of the aqueous 
eDNA particle size. 
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Figure 1. Diagram of particle size fractionation by sequential filtration. Surface water from 
throughout the water body was pumped into a 12 L carboy. 250 mL samples from this carboy 
were size fractionated by sequential filtration through filters of decreasing pore size, 180 - 0.2 
um. For each water body ten 250 mL samples provided filter retentate used in gravimetric 
analysis while ten samples were used for eDNA extraction and eDNA assays. Samples used for 
eDNA also provided a 15 mL subsample of 0.2 um filtrate from which eDNA was precipitated 
and pelleted by centrifugation. Only 150 mL of 1.0 um filtrate was used for filtration through the 
0.2 um filter to avoid filter clogging. 
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Figure 2. Percent of total carp eDNA, total eDNA, and SPM in each particle size fraction per 
water sample from the pond (A, B, C, respectively) and lake (D, E, F, respectively). *The final 
fractionation step was DNA precipitation using the 0.2 um filtrate. Filled symbols are winter; 
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open symbols are spring (only the lake was sampled in spring; SPM was not measured in spring). 
Data points are horizontally 'jittered' within each size fraction to reduce visual overlap of data 
points. Grey bars show the average value for each size fraction, which includes the supplemental 
spring data from the lake in (D) and (E). 
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Isocline Equation 



V, = 



V * E 



where. 



E ni = 0.788 

A = 11.85 
k = 0.351 



Figure 3. Isoclines showing combinations of filter pore size and water volume where the particle 
size distribution (PSD) of Carp eDNA predicts identical amounts of Carp eDNA captured. 
Isoclines are shown for five hypothetical examples of maximum throughput water volume for a 
0.2 um filter pore size (V0.2). The equation for calculating isoclines is shown using the 0.2 urn 
capture efficiency (E0.2) we estimated for Carp eDNA by fitting a Weibull CCDF to our 
cumulative size fractionation data; the scale (X) and shape (k) parameters from the Weibull 
CCDF model of Carp eDNA PSD are also used in this equation (see Results and Appendix SI 
for details). 
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Supporting Figure SI. Cumulative size fractionation data and the PSD models (Weibull CCDF) 
fitted to these data. Pooled lake and pond data are shown for Carp eDNA (A and E), total eDNA 
(B and F), and SPM (C and G). Panels (D) and (H) show the PSD models for each type of 
aquatic matter plotted together, without data points. All plots are shown with a logio x-axis (A, 
B, C, D) and a linear x-axis (E, F, G, H) to allow visualization of both the model fits and the PSD 
shapes. The estimated scale and shape parameters (k and k) from the Weibull CCDF model are 
shown on the plots. Supplemental spring data from the lake were not used for PSD modeling 
because SPM was not measured in spring. 
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Supporting Figure S2. Histograms of qPCR-measured carp eDNA copies per reaction. Note 



that the x-axis differs between pond and lake histograms. 
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Supporting Figure S3. Combinations of filter pore size and water volume used by previous 
researchers to capture aqueous eDNA from macrofauna. The one filled circle (20 urn pore size) 
represents a study where the filtered water volume was not reported. Note that filter material and 
surface area varied across studies. See Table SI for additional details and published sources. 
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Supporting Table SI. Published methods for capturing aqueous eDNA from macrofauna. 
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Supporting Appendix SI: Estimating model parameters for cumulative PSDs 

When modeling the cumulative particle size distribution (PSD) resulting from sequential 
filtration size fractionation, the explanatory variable (fixed and known) is the pore size of the 
filter. The response variable is the cumulative proportion of carp eDNA, total eDNA, or SPM 
retained at that pore size. Particle size fractions ranged from <0.2 um to > 180 um and the 
cumulative proportion necessarily ranges from zero to one. Because the final size fraction for 
each type of aquatic matter (<0.2 um for Carp eDNA and total eDNA, 0.2 - 1.0 um for SPM) 
necessarily had the cumulative proportion of one, these values were excluded from model 
estimation. 

Because there was no prior data regarding the size distribution of aqueous macrobial eDNA we 
chose to model it using a complementary cumulative distribution function (CCDF) that is 
flexible and has been widely used for modeling PSDs of other materials - the Weibull or 
Rosin/Rammler distribution. We used the CCDF of the Weibull distribution, 

(x\ k x > 0 

F(x; X,k) = e"w , where X > 0, 

k> 0 

To estimate the model parameter values using maximum likelihood with beta distributed errors 
(ML-P) we used custom code in Mathematica 8 (Wolfram Research, Champaign, IL, USA). We 
selected a beta error distribution because the random variable (proportion of total) is continuous 
and bound between zero and one (Ferrari and Cribari-Neto 2004): 
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0 < y < 1 

f(y> K & = W) r«i-M)0) y ^" 1(1 " y) (1 " M) * _1 . where 0 < M < 

0 > 0 

The gamma function is defined as T(z) = J 0 °° t z_1 e _t dt. 

From this model formulation we have a likelihood function of the form, 



1 I l-( e far) I ]</> 1 



Lafe,0;y,X) = nr=i- — V 7 (1-*)^ v 7/ , where 



r| | e~(i) )</> |r 



1-1 e (a) j U 

V \ \ // y 

« is the number observations. Previous research has used the beta distribution with logistic 
function for the mean (Polando et al. 2013). 

The study had twenty replicate size fractionations (10 pond, 10 lake) and the purpose of the 
model fitting was to find the best fit parameters under the twenty replicate samples. Thus the 
likelihood was fit using the multiplication of twenty likelihoods from independent samples 
(Edwards 1984), such that, 

L(A, k, 0; Y lt X lt Y 2 ,X 2 , ... Y 20 , X 2Q ) = L{X, k, 0; Y lt XJ * L(A, k, 0; Y 2 , X 2 ) * ... L{X, k, 0; Y 20 , X 20 ). 
All likelihood calculations were conducted in Mathematica 8 (Wolfram Research, Inc. 2010) and 
maximums were verified graphically. 
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